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Abstract 

The Casimir force is calculated analytically for configurations of two 
parallel plates and a spherical lens (sphere) above a plate with account 
of nonzero temperature, finite conductivity of the boundary metal and 
surface roughness. The permittivity of the metal is described by the plasma 
model. It is proved that in case of the plasma model the scattering 
formalism of quantum field theory in Matsubara formulation underlying 
Lifshitz formula is well defined and no modifications are needed concerning 
the zero-frequency contribution. The temperature correction to the Casimir 
force is found completely with respect to temperature and perturbatively (up 
to the second order in the relative penetration depth of electromagnetic zero- 
point oscillations into the metal) with respect to finite conductivity. The 
asymptotics of low and high temperatures are presented and contributions of 
longitudinal and perpendicular modes are determined separately. Serving as 
an example, aluminium test bodies are considered showing good agreement 
between the obtained analytical results and previously performed numerical 
computations. The roughness correction is formally included and formulas 
are given permitting to calculate the Casimir force under the influence of all 
relevant factors. 
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1 Introduction 

Lately considerable progress had been made both in theoretical and 
experimental investigation of the Casimir effect. This effect predicted by 
H.B.G. Casimir 1 more than fifty years ago, consists in the interaction of two 
neutral, conducting bodies placed in vacuum close to each other. The Casimir 
effect results from the disturbance by the conducting boundaries of the zero- 
point electromagnetic oscillations. It plays an important role in various fields 
of physics such as elementary particle theory, condensed matter physics, 
atomic physics, gravitation and cosmology, and stimulated new investigations 
in mathematical physics (see the monographs 2-5 ). Recently the Casimir effect 
found applications 6 " 11 for obtaining rather strong constraints on hypothetical 
long-range interactions inspired by the physics of extra dimensions, by unified 
gauge theories, supersymmetry and supergravity. Furthermore, topical 
nanoelectromechanical devices were proposed 12 " 14 which are based on the 
use of the Casimir force. 

In precision experiments on the measurement of the Casimir force 14-19 
different influential factors must be accounted for, such as nonzero 
temperature, finite conductivity of the boundary metal and surface 
roughness. Theoretically, each factor was investigated in a number of papers 
(see, e.g., Refs. 20-24 for the influence of nonzero temperature, Refs. 22, 25-29 
for the role of finite conductivity and Refs. 30-36 for the surface roughness). 
The combined effect of different corrections was discussed in Refs. 23, 24, 35 
(for a detailed discussion of this subject see the recent review 37 ). 

Investigation of the combined effect of nonzero temperature along with the 
finite conductivity of the boundary metal proved to be the most complicated 
task leading to controversial results. The starting point to theoretically 
describe this effect is the Lifshitz theory 38 originally being developed for 
dielectrics. In order to describe on the base of Lifshitz' theory the Casimir 
force between plates made of an ideal metal a special prescription was 



4 



suggested in Ref. 22 which demands to consider first the limit of infinite 
dielectric permittivity before setting the frequency equal to zero. Using this 
prescription, the results of the Lifshitz theory for ideal metals agreed with 
the results obtained by the application of quantum field theory using the 
idealized boundary conditions. 20 ' 21 

In the last year, several authors attempted to apply the Lifshitz theory to 
calculate the Casimir force between plates made of real metals. In Ref. 23 the 
plasma model was used to describe the dependence of dielectric permittivity 
on frequency. In Ref. 39 the Drude dielectric function (being a generalization 
of the plasma model by taking into account the relaxation processes) was 
substituted into the Lifshitz formula. However, in the limit of zero relaxation 
frequency the results of Ref. 39 do not coincide with those of Ref. 23 although 
Drude's model turns into plasma model in this limit. In Ref. 24 the results of 
Ref. 23 obtained using the plasma model were independently confirmed and 
doubts were casted on the calculations of Ref. 39 using the Drude model. 
It was noticed 24 that the high temperature asymptotics of the Casimir force 
between real metals computed in Ref. 39 is two times smaller than in the 
case of an ideal metal - - independently of how high the conductivity of the 
real metal is - - which is a nonphysical property. As was noted in Ref. 40, 
the computations of Ref. 39 are also in contradiction with the experiment. 15 

To improve this situation in Ref. 41 a detailed computation of the Casimir 
force at nonzero temperature were performed based on Drude's model and a 
Lifshitz formula with some modified zero-frequency term. The modification 
made is based on a generalization of the prescription of Ref. 22 for the case 
of real metals. The results of Ref. 41 are in agreement with both limiting 
cases of an ideal metal and a metal described by the plasma model. In 
Refs. 42,43 one more result was obtained for the temperature Casimir force 
between real metals which, however, disagrees with both the results of Ref. 39 
from one side and of Refs. 23, 24 from the other. According to Refs. 42, 43 at 
small frequencies all real metals are indistinguishable from the ideal metal. 
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This leads to the absence of any finite conductivity correction to the Casimir 
force starting from moderate separations of several micrometers between the 
test bodies. Even more, in the approach of Refs. 42,43 this property is 
independent on the quality of the real metal. 

As is seen from the above, at the time being there is no agreement in the 
literature concerning the calculation of the Casimir force acting between real 
metals at nonzero temperature. Different results are obtained on this subject 
by different authors starting from one and the same theoretical foundation 
given by the Lifshitz formula for dielectrics. In the present paper we discuss 
the scattering formalism of quantum field theory at finite temperature in the 
Matsubara formulation underlying Lifshitz formula and we argue that it leads 
to well defined and consistent results both physically and mathematically 
if the dielectric function is described by the plasma model. This gives 
additional support to the results of Refs. 23, 24. Contrary, when the dielectric 
function of the Drude model is used the scattering formalism becomes 
inconsistent causing the nonphysical results of Refs. 39,42,43. In this case 
the modification of Lifshitz formula suggested in Ref. 41 is needed. 

Below, starting from the dielectric function of the plasma model we 
calculate the Casimir force for the configurations of two parallel plates 
and a spherical lens (sphere) above a plate made of real metals. The 
temperature corrections are taken into account completely and the effects 
of finite conductivity of the boundary metal are treated perturbatively up 
to second order in some small parameter having the meaning of the relative 
penetration depth of electromagnetic zero-point oscillations into the metal. 
Note that in the analytical computations of the previous paper 23 the finite 
conductivity corrections were calculated up to the first order only, whereas in 
Ref. 24 both corrections (due to nonzero temperature and finite conductivity) 
were treated perturbatively. The analytical results obtained below are 
compared with the results of numerical computations and good agreement 
is observed for all space separations exceeding the plasma wavelength of the 
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boundary metal. We also include roughness correction and demonstrate a 
way how to take into account the influence of real conditions (which include 
all three types of corrections) onto the Casimir force. 

The paper is organized as follows. In Sec. 2 the general formalism is 
presented and the scopes of its consistency are elucidated. In Sec. 3 the 
temperature correction to the Casimir force is calculated for the configuration 
of two parallel plates up to the second perturbation order in the relative 
penetration depth. Sec. 4 contains the analogous results for the configuration 
of a sphere (spherical lens) above a plate. In Sec. 5 the roughness correction 
is taken into account along with nonzero temperature and finite conductivity 
corrections. In Sec. 6 the reader finds our conclusions and discussion. 

2 General formalism and its scopes 

Lifshitz' original derivation 38 of the Casimir force at nonzero temperature 
acting between two dielectric semispaces separated by a gap was based 
on the assumption that the dielectric materials can be considered as 
continuous media characterized by randomly fluctuating sources. The 
modern derivation 37 is based on quantum field theory at nonzero temperature, 
T ^ 0, in the Matsubara formulation. Thereby one considers the Euclidean 
version with the electromagnetic field periodic in the Euclidean time variable 
within the time interval (3 = 7&/(&bT), where Ub and Ti are the Boltzmann 
and Planck constants, respectively. 

Let us consider two dielectric semispaces with frequency dependent 
permittivity e{uS) restricted by two planes at z = ±a/2 and separated 
by a vacuum gap of width a between them. Due to periodicity in the 
time coordinate, the frequency spectrum is discrete uj\ = 27r///3, whith 
/ = ...— 2, —1, 0, 1, 2, ... . The calculation of the free energy is reduced 
to the solution of a one-dimensional scattering problem on the z-axis. In fact 
an electromagnetic wave which is coming from the left, or from the right, in 
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the dielectric semispaces will be scattered on the vacuum gap and there will 
be a reflected and a transmitted wave. The free energy of the field per unit 
area, E 88 , is calculated with the help of (-regularization method. The result 



is 
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where s xl fejj and fej_) are the scattering coefficients for parallel 

and perpendicular polarizations, respectively, k± = (k x , k y ) is the wave vector 
in the planes perpendicular to the z-axis, and = 2irl/(3. The solution of 
the scattering problem reads 37 

H 4e(i6)fci«e*» a 



'ii 



'11 



qia 



4kiqie 



ha 



(qi + k) - ( qi - k) e-u a 
where the following notations are introduced 



(2) 



% + k 2 



\ c 2 



\ 



£2 



(3) 



Qi = 

\ cr 

Now we substitute (@) into ([I]) and perform renormalization in order to get 
the free energy equal to zero in the case of infinitely far remote plates. 28 In 
terms of reflection coefficients rn, r± for the electromagnetic waves of the two 
different polarizations one obtains 
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The Casimir force per unit area acting between two semispaces is obtained 
as —dE ss /da with the result 



k B T °r 
F ss (a) = — - — / k ± dk ± qi 
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(6) 



This equation, up to change of variables, coincides with the original Lifshitz 
expression 38 for the Casimir force between dielectrics at nonzero temperature. 

Our aim is to apply Eqs. (|])-(@) in the case of test bodies made of 
real metals. Then, the zero-frequency contribution for £o = in Eq. (||) 
may become indefinite in the case of perpendicular polarization when the 
dielectric permittivity turns into infinity. For example, let us consider the 
Drude dielectric function on the imaginary axis 



where cu p is the plasma frequency, 7 is the relaxation frequency, which 
gives a good approximation of dielectric properties for some metals, e.g., 
for aluminium. Note that for dielectrics — » £0 when £ — ► 0. It is evident 
that 



both for Drude's model and dielectrics. By this reason, it follows from Eq. (j3|) 
that qo = ko = k±. Strictly speaking, the mathematical derivation leading to 
Eq. (@) for sj\ is inapplicable in that case. Instead, the direct solution of the 
scattering problem in the case qo = ko gives the result that sj\ is arbitrary 
and — 0, where Su is the nondiagonal element of the scattering matrix. 
In the case of dielectrics the unitarity condition is valid which immediately 
leads to [s^] 2 = 1 and, due to dispersion relation, to sjj = 1. In fact, the 
same result is obtained from Eq. (0) in the limit qo = ko- However, as to 
the case of the Drude model, which describes a medium with dissipation, the 
unitarity condition is not applicable, and, therefore, the scattering coefficient 
sj\ remains indefinite. Because of this, the direct application of Lifshitz 
formula (as in Ref. 39) leads to incorrect results. To apply the Lifshitz formula 
at nonzero temperature in combination with Drude model some additional 
prescription is needed to give the definite value to the zero-frequency term 
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(see Ref. 41 for details). 

At the same time for the longitudinal polarization the scattering coefficient 
is well defined in the limit of zero frequency. In the case of the Drude model, 
up to terms independent of a, it has the limiting value 



'ii 



1 



-2ak i 



when £ — > 0, 
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i.e., the same as for ideal metals. For dielectrics 
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when £ — > 0. 
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Thus, the transition from Eq. (fj) to Eqs. (f|), (|6|) is unjustified in the case of 
Drude's dielectric function. This explains why the nonphysical results were 
obtained when substituting Eq. (0) into Eq. (§) (see Introduction). 
Another model of the dielectric function is the plasma one, 

e(<f) = 1 + (11) 



which is the limiting case of (0) when the relaxation frequency 7 goes to zero. 
In the case of the plasma dielectric function 



{ 2 e(i{)^^0 whenf-0. 



(12) 



As a consequence, here qo ^ ko and the limiting value of the perpendicular 
scattering coefficient is given by 

-1 



—2ak 1 



k 



_L 



when £ — > 0. 



(13) 



The limiting value of the longitudinal scattering coefficient in the case of 
plasma model is the same as in Eq. (g). Because of this, the scattering 
problem is well defined for the dielectric function ([TT|) and Eqs. ([|) and (§) 
can be reliably applied to calculate the free energy and the Casimir force. No 
additional prescriptions or modifications are admissible in the case of plasma 
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model. Because of this, the manipulations of Refs. 42,43 changing the zero- 
frequency term of Lifshitz formula for both plasma and Drude models seem 
to be unfounded. In the case of the plasma model they lead to disagreement 
with the results of Refs. 23, 24 where no modifications of Lifshitz formula 
have been made. 

In the next section the analytical computations of the temperature 
correction to the Casimir force are performed on the basis of Eqs. ([|), (|5|) by 
using the plasma model ([TT]) . 



3 Temperature correction to the Casimir force for two 
parallel plates made of real metal 



We start with the Lifshitz formula 
variables 

y = 2aqi = 2a 

resulting in 

k B T 



and rewrite it in terms of dimensionless 
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Here the reflection coefficients in new variables take the form 
ey - 
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ey + ^(e- l)xf + y 2 
and e = e [icxi/(2a)}. 

To separate within ( |T5| ) the contribution of temperature T = and the 
temperature correction one can use the representation of this formula in terms 
of continuous x instead of discrete summation in x\. Applying the Poisson 
summation formula one obtains from (|T5|) 20 ' 22 ' 23 
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where t = T e ff/T, k,BT e fj — he/ (2a). The reflection coefficients preserve 
their form with a change xi — >• x. 

In Eq. (|T7|) the term with / = is the Casimir force at zero temperature, 22 
whereas the terms with / ^ represent the temperature corrections. For 
further needs it is convenient to write the temperature correction of the 
Casimir force acting between real metals as a sum of longitudinal and 
perpendicular contributions, 



A T F ss (a) = A^F ss (a) + A^F ss (a), 



(18) 
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To compute the temperature correction ([ISO, (HI) one should use some 
specific functional dependence of the dielectric permittivity on the frequency. 
Here we use the plasma model (|TT|) for which the theory is well defined (see 
Sec. 2). In terms of dimensionless variables the dielectric function ( |TT|) is given 
by 



e = e(x) = 1 + 



1 + 



4a z 



(20) 



x^ SqX 2 ' 

where u p = 2acu p /c, 5o = c/u p is the effective penetration depth of 
electromagnetic zero-point oscillations into the real metal. For the space 
separations between the plates a ^> 5o the natural small parameter is 
^o/o < 1. In fact this condition is valid for a > \ p , where X p = 2ttc/u p 
is the effective plasma wavelength, since $o = X p /(2tt). 

Here we calculate the temperature correction QIBQ , fligQ analytically taking 
completely into account the nonzero temperature and using the perturbation 
theory up to the second power in the small parameter So/ a in order to take 
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approximate account of the finite conductivity of the boundary metal. To 
perform the computations let us expand first the expressions in Eq. 
containing reflection coefficients, in powers of 5$/ a. The result is 
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Substituting ( |2T| ) into and calculating integrals with the help of the 
formulas 3.951(12, 13) of Ref. 44 one obtains finally the contribution of the 
longitudinal mode 
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Quite analogically the contribution of the perpendicular mode is 
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Finally, the total temperature correction ([T8|) is given by 

1 7T 3 coth(7r/t) 
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Note that the first two terms on the right-hand side of ( |24| ) which are of zeroth 
order in 5o/a coincide with the well-known result for the ideal metal 21 ' 22 , 
whereas the coefficient of the first power in So/ a was first obtained in Ref. 23. 

In the case of low temperatures (small separations), t ^> 1, one can 
substitute the hyperbolic functions by their asymptotics. Preserving the 
largest of the exponentially small contributions and performing summations 
of the power ones we obtain from 
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where ((z) is Riemann's zeta function. It is seen that the second perturbative 
order in 5q/cl is exponentially small in t and does not contain purely powers 
in t contributions like the first order term of (|Z5p . This is in agreement with 
the perturbation results of Ref. 24 where the double perturbation theory 
in the small parameters So/ a and \/t was developed. At the same time, 
under the natural supposition (So/a)t ~ 1 the second order term turns out to 
be approximately 90 times larger than the exponentially small contribution 
in the zeroth order term and 7 times larger than the exponentially small 
contribution in the first order term. 
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Now consider the case of high temperatures (large separations) when t <C 
1. It is more simple to extract it not from Eq. ([24] ) but directly from Eqs. fllSQ, 
(|T91) , (]2Tj) . To do this one should perform the integration with respect to x 
in the same way as above and then change the order of summation and 
integrations with respect to y. Due to the smallness of t all summations can 
be performed by the use of the formula 44 

sin(lty) 7r — ty 



E 



I 



(26) 



which is valid for < ty < 2n. In the further integrations with respect to 
y all functions under the integrals decrease with y as exp(— y), so that the 
infinite upper limit of the integration can be changed for y = {2ir/t) — a, 
where a > 0, with the required accuracy. As a result, the high temperature 
limit of the temperature correction to the Casimir force between real metals 
is given by 
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Let us discuss the application range of the analytical result and the 
asymptotic representations ( |25| ) and (|27|). Bearing in mind that according to 
Eq. ( \T7\ ) the Casimir force at nonzero temperature and finite conductivity is 
given by 

F ss (a) = F ss (a; T = 0) + A T F ss (a), (28) 
it is convenient to compute the quantity 

" F„(a;T = 0) V ' 

Then the value of (1 + k ss ) has the meaning of a correction factor. Indeed, 
multiplying the Casimir force F ss (a; T = 0) computed with account of finite 
conductivity at zero temperature by (1 + k ss ) one obtains the Casimir force 
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at both nonzero temperature and finite conductivity. Note that F ss (a; T = 0) 
was computed in a number of papers 27-29 , and below we use the numerical 
and analytical results obtained there to calculate k ss . 

In Table 1 the values of k ss are presented at several separations computed 
by the use of Eq. (|Z4Q (second column), by the use of low-temperature 
asymptotics (third column), and high-temperature asymptotics (fourth 
column). In all computations, as an example, the value of the plasma 
frequency tu p = 12.5 eV is used as for aluminium 45 (this corresponds to the 
plasma wavelength of approximately X p ~ 99 nm). It is seen from Table 1 that 
the asymptotic of low temperatures gives the same values of k ss as Eq. (|Z3Q 
at all separations a < 2 /xm (for a > 3/im it is not applicable). Comparing 
data of columns two and four one can conclude that the asymptotics of high 
temperatures works good for a > 7 /im and is not applicable for a < 5 /im. In 
the transition region 3 /im < a < 5 /im neither of the asymptotics but Eq. (|2j|) 
itself should be used to compute the temperature correction to the Casimir 
force acting between real metals. It is noticeable, that data of column 2 are 
practically the same irrespective of whether one or two perturbation orders 
are taken into account. These data coincide also with the results of numerical 
computations by Eqs. flTBl)-( EZUl) in all separation range 0.1 /im < a < 10 /im. 

It is interesting also to discuss the comparative contribution to the 
temperature correction which results from the longitudinal and perpendicular 
modes given by Eqs. (|22|), (P3|). At small separations the contribution of 
A^F SS to the temperature correction dominates the contribution of Aj^F ss . 
By way of example, at a = 0.1 fim one has A^F SS / A^F SS = 43.1. This ratio, 
however, quickly decreases with the increase of space separation. Thus, at 
a = 0.3 /im it is equal to 5.67, whereas at a = 0.5 /im and a = 0.7 /im it 
equals to 2.68 and 1.86, respectively. At large separations both modes give 
almost equivalent contribution to the temperature correction. For example, 
at a = 7 /mi the abovementioned ratio is equal to 1.015 and at a = 10 /im it 
is equal to 1.01. 
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4 Temperature correction for a sphere above a plate 
made of real metal 

In most of the experiments on measurement of the Casimir force the 
configuration of a sphere (spherical lens) placed above a plate (semispace) is 
used 15-18 (in fact the configuration of two crossed cylinders, as in Ref. 19, is 
equivalent to it). The expression for the Casimir force at nonzero temperature 
acting in this configuration can be obtained by means of the proximity force 
theorem 46 

F sl (a) = 27rRE ss (a), (30) 

where E ss (a) is the free energy per unit area of the two plates defined in 
Eq. (f|), R is the curvature radius of the sphere (spherical lens). In terms of 
the dimensionless variables of Eq. ([T|]) the force acting between a lens and a 
plate is 

F 8 i(a) = E / ydy {in [l - rjj (x h y) e~ y 

+ \n[l-rl( Xl ,y)e-y]}, (31) 

where the reflection coefficients run) are defined in ([Hf). 

After applying the Poisson summation formula, Eq. (£JT]) can be rewritten 
in the form analogical to Eq. ([17]) 

F s i(a) = J2 J ydy J dxcos(ltx) 

x {In [1 - r\ (x, y) e~ y ] + In [l - r| (x, y) e~ y ] } . (32) 

Once more, the term with / = is the Casimir force at zero temperature, the 
terms with / ^ represent the temperature corrections to it. In accordance 
with Eq. ( |I8|) the temperature correction can be splitted into a sum of 
longitudinal and perpendicular contributions (with a change of index ss — > si) 
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expressed by 
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To compute the temperature correction analytically by using the plasma 
model we substitute Eq. ( |2D| ) into Eq. ( |33Q and, again, expand into powers 
of the small parameter 8o/a 



In 1 — r 2 | (x, y) e v 
- 2 



On 



a 



ln(l 



+ 2- 



3" 



a y(e y — 1) 



y\ev - I 



+ 



'A3' 

a 3 



(34) 



In 1 — r\ (x, y) e 



ln(l 



+ 2 



<5q 2/ 
a ev - 1 



2 - 

V a 



?/ 2 e y 



r 



+ o,f 



Let us calculate first the longitudinal temperature correction which is 
obtained by the substitution of the first equality from fl3jj ) into (|33|) . All 
integrals with respect to x are trivial. The resulting integrals with respect 
to y can be found in Ref. 44 (formulas 3.951(12, 13)) except of the following 
one which, however, also can be computed analytically: 
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coth(7r/£) — 

4 V ' 47T 



7T 



h 27rZtln(l - e 

6 v 



-2nlt 



2ttW 

2ttU _ 1 



Li 2 (e 



-2wlt 



where Li 2 (^) is the polylogarithm function. As a result the contribution of 
the longitudinal modes to the temperature correction is 



AJ^,,(a) = 



hcR ™ 



7T 



87ra 3 rri[2(/t) 3 



+ 2- 



5c 



a 



71 



(ity 



coth(7r/t) 



coth(7r/t) 

3 IT 2 
+ 



1 



+ 



7T 



1 



(lt) A 2{lt) 2 smh 2 {7rlt) 
7T 3 coth(7r/t) 



(36) 



+ 



[Ity {lt) 2 smh 2 (irlt) ltsmh 2 (irlt) 
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a 



7T 



7T 



6 



7r(lty 



(It) 5 smh 2 (7rlt) 
27r^ln(l - e- 2nlt 



1 - 3coth 2 (7r/f) 

27T 2 (lt) 2 



coth(7r/t) 



,2nlt 



nit 

Li 2 (e- 2 * lt ) 



e-"" - 1 

The contribution of the perpendicular modes is calculated simply as 

^ 2 1 



. i , , hcR ~ 



7T 



[2(/f) ; 



coth(7r/t) 



1 7T 

+ 



{Ity 2{lt) 2 smh 2 (nlt) 



+ 2^ 
a 



7T 3 COth(TT^) 1 

¥sinh 2 (7r/t) ~ (/^ 



(37) 



7T 



[3 coth(Tr^) + ?rZt - 3tt/£ coth 2 (7r/£)] 



U/ sinh 2 (7r/t) 

Putting together the contributions of both modes the total temperature 
correction for the configuration of a lens above a plate is obtained: 



%cR 00 [ 7T 



1 



+ 



7T 



1 



a 



IT 



(ity 



coth(7r/t) 



+ 



71" 



(Ity 2{lt) 2 sinh 2 {nlt) 
2tt 3 coth(7r/t) 



(38) 



1 



+ 



'So\ 2 
, a J 

6 

^njti)' 



71 



+ 



(/£) 4 ' (/f) 2 sinh 2 (7r/f) It sinh 2 (nlt) 
2tt 4 / . . o , . coth(irlt) 



(It) 5 smh 2 {nlt) 



2nltln 1- e 



1 - 3coth 2 (7r/f) + 
2n 2 {lt) 2 



1 



Tilt 

Li 2 (e" 2 *" 



(?rZt) : 



e*"" - 1 

The terms of zeroth order in Sq/ 'a in the right-hand side of ([38|) coincide with 
the known result for an ideal metal. 37 The coefficient of the first power in 
5o/a was already obtained in Ref. 23. 

Now, let us consider the limiting cases of Eq. (|38|) corresponding to low 
and high temperatures (small and large separations). At low temperatures, 
t ^> 1, and, preserving the largest of the exponentially small contributions, 
one obtains from Eq. (|35| ) 



A T F 8l (a) 



TicR f ttC(3) tt 4 



47ra 3 



+- 



Sr. 



a 



C(3) 



2£ 3 
2tt 4 
45£ 4 



+ 



90£ 4 
8tt 3 



+ 



2^ 

I 2 



-2nt 



(39) 



-2nt 



t 



%\ 2 



a 



<(5) 
t 5 



167rV 27rt 
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It is noticeable, that for the configuration of a lens above a plate the second 
perturbative order in Sq/cl contains power-type contributions in t, not only 
exponentially small ones. It is, however, of order t~ 5 in agreement with 
Ref. 24 where the absence of temperature corrections of powers lower than 
1/t 5 was proved for the perturbation orders (5o/a) k with k = 2, 3, 4, 5, 6. 

The limit of high temperatures can be obtained in the same way as in 
Sec. 3, i.e., starting from Eq. fl33Q and using Eq. (|2^). The result is 

A T F sl (a) « I 40 

v ; 47ra 3 I 2t 90 y J 



a 



<(3) 2tt 



2 

+ 



3ttC(3) 4tt 



t 25 



t 45 

To determine the range of applicability of both asymptotic representations 
we compute the quantity 

A T F d (a) 

k " = F sl (a;T = 0y (41) 
where F s i(a; T = 0) is the force at zero temperature calculated with account of 
finite conductivity. 27-29 The value of (l + fe s /) has the meaning of a correction 
factor to it. The total Casimir force acting between real metals at nonzero 
temperature is given by 

F sl (a) = (l + k al )F al (a;T = 0). (42) 

In Table 2 the values of k s i are presented for aluminium computed (i) 
by Eq. (|38|) (second column), (ii) by the low-temperature asymptotic (third 
column), and (iii) by the high-temperature asymptotic (fourth column). In 
analogy with the case of two parallel plates, the asymptotics of low and 
high temperatures work good at separations a < 2 jum and a > 6 /mi, 
respectively. The data of column 2 are practically the same in first- and 
second-order perturbation theory. They coincide also with the results of 
numerical computations by the use of Eq. fl33f) . 

The comparative contribution of the longitudinal and perpendicular modes 
to the temperature correction to the Casimir force (see Eqs. (|35Q , (|37|) ) is 
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different in comparision with the case of two plates. Here, A^F s i/ Aj^F s i = 
1.64 at a = 0.1 /mi, 1.21 at a = 0.3 /im, and 1.13 at a = 0.5 /im. This ratio 
decreases slowly to the value 1.006 at a = 10 /im. Thus, the contributions of 
both modes are approximately equal to each other at all separations (for two 
plates the contribution of the longitudinal mode significantly dominates at 
smallest separations). This is explained by the presence of the term ~ in 
the zeroth order contribution to AtF s i (in the case of two plates the zeroth 
order contribution is ~ t~ A ). 



5 Casimir force with account of roughness 



The above Eqs. fl28|), flfflp give us the Casimir force computed at nonzero 



temperature with account of finite conductivity of the boundary metal. 
Except for the finite value of conductivity, real metallic boundaries are 
characterized also by some surface roughness. According to the results of 
Ref. 35 obtained at zero temperature, for a wide range of surface roughness 
it can be taken into account by some kind of geometrical approach using the 
averaging of the Casimir force over the rough surface. Here, we generalize 
this approach for the case of nonzero temperature. 

Let two large metallic plates of dimension L x L be covered by small 
roughness. Then the distance between two points of the boundary surfaces 
of different plates with the coordinates (x, y) can be expressed as 

a(x,y) = a + f(x,y). (43) 

Here, f(x, y) is simply expressed by the functions describing roughness on 
both surfaces. The mean distance between the plates ao is defined in such a 
way that 

L L 

J dx J dyf{x,y) = 0. (44) 

—L -L 

As a result, the Casimir force taking all real conditions into account (i.e. 
nonzero temperature, finite conductivity of a metal and surface roughness) 
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is given as 




Remind that F ss here is given by Eq. flZ8Q with a change of separation distance 
a for the one defined in Eq. (||3|). The analogical result can be obtained for 
the configuration of a sphere above a plate by the use of the proximity force 
theorem. 46 As is shown in Ref. 37, Eq. fl35j) gives the same results as the more 
fundamental methods for accounting roughness, e.g., based on the specific 
forms of interatomic potentials or Green's function method. 

Finite conductivity and temperature corrections to the Casimir force were 
computed in Sees. 3, 4 in the whole distance range 0.1 /im < a < 10 /mi. 
Surface roughness makes the most important contribution to the Casimir 
force for separation distances a < 1 /im. For such distances (and also for 
1 /im < a < 3 /im) the asymptotic of low temperatures of Eq. (|Z5p can be 
substituted into Eq. ([|5]) to calculate the Casimir force under real conditions. 
As to the transition region 3 /im < a < 5 /im, the more exact Eq. (|24] ) should 
be used there (at larger separations roughness corrections are negligible). 

6 Conclusions and discussion 

As was argued above, the dielectric function as it results from the plasma 
model can be reliably used to calculate the Casimir force acting between real 
metals at nonzero temperature. The scattering theory underlying the Lifshitz 
formula is well defined in the case of the plasma model and its application is 
straightforward. No additional prescriptions are needed like those formulated 
in Ref. 22 for ideal metal or in Ref. 41 for real metals described by the Drude 
dielectric function. 

We calculated the temperature correction to the Casimir force between 
real metals in the configuration of two parallel plates (two semispaces) and 
for a lens (sphere) above a plate. The analytical expressions for these 
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corrections were obtained which are exact with respect to the temperature 
but perturbative with respect to the effects of finite conductivity. These 
effects were taken into account up to the second order in a small parameter 
having the meaning of the relative penetration depth of electromagnetic 
zero-point oscillations into the metal. The asymptotics of the obtained 
expressions were presented at both low and high temperatures relative to 
kBT e ff = he/ (2a). The asymptotical formulas are in good agreement with 
the exact ones except of a narrow transition region between the cases of low 
and high temperatures. The scopes of these regions are determined. In a wide 
separation range from 0.1 fim till 10 fim the obtained analytical results are 
in perfect accordance with the results of numerical computations performed 
earlier. The comparative contributions of the longitudinal and perpendicular 
modes to the temperature correction were determined. Modification of the 
obtained results taking the surface roughness into account was given. This 
together permits to evaluate the Casimir force under the influence of real 
conditions which include nonzero temperature, finite conductivity of the 
boundary metal and surface roughness with a precision of several percent. 

The above results are topical ones for the precision measurements of the 
Casimir force. In view of fundamental and technological applications of the 
Casimir effect mentioned in Introduction, there is a great theoretical challenge 
to account for real experimental conditions. In fact, a theory is required which 
makes it possible to calculate the Casimir force with a precision of being 
better than one percent. For this purpose one should especially examine the 
optical properties of the test bodies in use, their surface roughness and take 
into account the spatial dispersion (in the case when there are thin layers 
covering the test bodies). The more precise future theory should take into 
account also the effects of dissipation which are neglected in the presently 
used Lifshitz formula. 
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Table 1: The values of k ss computed by the use of Eq. ( |2"4] ) in comparison with the asymptotic 
values at low [Eq. and high [Eq. (^)] temperatures. 



a k ss computed k ss at low k ss at high 
(/im) by Eq. ( ^4|) temperatures temperatures 



0.1 


6.58x10- 


-G 


6.58x10- 


-G 




0.3 


5.48x10" 


-5 


5.48x10- 


-5 




0.5 


2.12x10" 


-4 


2.12x10- 


-4 




0.7 


6.05x10- 


-4 


6.05x10- 


-4 




0.9 


1.42x10- 


-3 


1.42x10- 


-3 




2.0 


2.74x10- 


-2 


2.76x10- 


-2 




3.0 


0.1228 




0.1339 






4.0 


0.3100 




0.4159 




0.2302 


5.0 


0.5630 








0.5349 


6.0 


0.8487 








0.8400 


7.0 


1.147 








1.144 


8.0 


1.450 








1.449 


10 


2.059 








2.059 
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Table 2: The values of k s i computed by the use of Eq. in comparison with the asymptotic 
values at low [Eq. and high [Eq. (fEJ)] temperatures. 



k s i computed k s i at low 



k s i at high 



(/im) by Eq. (j38|) temperatures temperatures 



0.1 


6.69x10 


-5 


6.70x10- 


-5 




0.3 


1.08x10" 


-3 


1.08x10- 


-3 




0.5 


4.33x10- 


-3 


4.33x10- 


-3 




0.7 


1.09x10" 


-2 


1.09x10- 


-2 




0.9 


2.181x10 


-2 


2.181x10 


-2 




2.0 


0.1829 




0.1828 






3.0 


0.4810 




0.4762 




0.3857 


4.0 


0.8736 




0.8113 




0.8428 


5.0 


1.309 








1.300 


6.0 


1.759 








1.757 


7.0 


2.215 








2.214 


8.0 


2.671 








2.671 


10 


3.585 








3.585 



